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ABSTRACT 


The  initial  resultsof  a  variable  amplitude  fatigue  crack  growth 
rate  test  program  are  presented.  The  test  program  was  designed  to 
study  the  Influence  that  the  stress  intensity  factor  gradient  (dK/da) 
Had'Bh  crack  growth  rate  behavior  in  7075-T6  aluminum  aTr(^^  Constant 
stress  intensity  factor- conditions  were  used  to  study  this  behavior, 
using  a  twenty-four  inch  wide  center-cracked  panel. 

Based  on  earlier  work  by  Hillberry  on  constant  amplitude  ,crae]^  ^ 
growth  rate  behavior,  a  three-parameter  lognormal  distribution  function 
was  selected  to  describe  the  crack  growth  rate  data.  Comparison  of 
crack  lengths  showed  no  gradient  effect.  Parameters  of  the  statistical 
distribution  were  compared  using  the  secant  and  incremental  polynomial 
methods  of  differentiation  and  for  different  flight  intervals  achieved 
by  removing  points  prior  to  the  differentiation  process. 
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SECTION  I 


INTRODUCTION 

Fllght-by-f light  crack  growth  rate  (Aa/AF)  data  are  important  to 
the  structural  design  of  fatigue  critical  aerospace  components.  This 
paper  documents  an  investigation  of  the  variability  associated  with 
flight-by-flight  crack  growth  rate  behavior.  Variability  in  crack 
growth  rate  data  occurs  because  of  error  in  crack  length  measurement, 
lack  of  precise  load  and  environmental  control,  as  well  as  variability 
in  material  behavior.  The  type  of  differentiation  technique  employed 
also  effects  the  reported  variability.  Measurement,  load,  and  environ¬ 
mental  discrepancies  can  be  minimized  by  standardizing  test  procedures. 

For  example,  choice  of  the  proper  crack  length  measurement  interval 
will  help  to  minimize  the  relative  size  of  the  measurement  error  and 
thereby  decrease  the  reported  variability.  The  influence  of  variations 
in  test  parameters  should  be  realized  when  employing  crack  growth  rate 
test  data  for  design  purposes,  as  reported  by  Hillberry  ,  Wei  and 
others. 

In  order  to  establish  baseline  crack  growth  rate  data  for  this 
study,  a  variable  amplitude  load  history  was  applied  to  a  twenty-four 
inch  center-cracked  panel.  From  earlier  experience,  this  particular 
loading  sequence  was  known  to  produce  steady  state  crack  growth  rate 
behavior^^’^\  Steady  state  behavior  is  similar  to  that  experienced 
under  constant  amplitude  fatigue  loading. 

A  "short"  load  history  was  chosen  so  that  the  influence  of  the 
stress  intensity  factor  gradient  (dK/da)  on  the  crack  growth  rate  behavior 


in  7075— T6  aluminum  alloy  could  be  studied.  Load  shedding  techniques 
were  employed  to  obtain  constant  stress  intensity  factor  conditions  to 
study  the  resulting  variability.  A  "short"  load  history  has  been  defined 
by  Elber^^^  as  one  in  which  the  crack  growth  during  one  repeat  load 
interval  is  less  than  the  plastic  zone  size. 

In  this  investigation,  flight-by-flight  crack  growth  rate  data 
resulting  from  the  secant  method  of  differentiation  were  fit  to  four 
distribution  functions.  This  was  done  to  select  the  function  which  best 
describes  the  distribution  of  the  Aa/AF  data.  A  measure  of  the  variability 
was  then  obtained  in  terms  of  the  distribution  parameters.  It  was  first 
necessary  to  show  .that  the  variability  did  not  change  as  the  crack 
advanced  under  constant  stress  intensity  factor  conditions.  Increasing 
the  crack  length  measurement  interval  substantially  decreased  the 
variability. 
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SECTION  2 


TEST  DESCRIPTION 


2.1  Load  History 

A  modified  design  stress  history  for  bomber  type  aircraft, 
containing  135,000  cycles  per  lifetime,  was  used  for  this  test  program. 

The  original  stress  history  was  composed  of  three  mission  blocks  as 
studied  by  Potter The  first  contained  only  once  per  flight  load 
events.  The  second  and  third  contained  all  load  events  of  the  first 
plus  additional  load  events  which  occurred  once  in  ten  flights  and 
once  in  one  hundred  flights,  respectively.  For  the  purpose  of  establish¬ 
ing  steady-state  crack  growth  rate  behavior,  the  higher  periodic 
loads  were  omitted;  only  the  once  per  flight  load  events  were  included 
in  the  initial  repetitive  sequence.  All  compression  loads  were  taken 
out  of  the  load  history  (see  Figure  1) . 

The  load  history  was  such  that  the  crack  never  grew  out  of  the 
plastic  zone  in  a  single  flight.  The  highest  load  appearing  in  the 
modified  load  history  was  88  percent  of  the  highest  load  level  v/hich  occurred 
in  the  original  history.  Because  of  the  unique  relation  between  loads 
in  the  load  history,  constant  stress  intensity  factor  conditions  were 
accomplished  through  the  reduction  of  the  maximum  repeating  load  as 
suggested  by  the  Fedderson  secant  formula^^^ .  This  insured  that  the 
maximum  stress  intensity  factor,  K  ,  was  kept  constant  within  one 

TlldX 

percent  as  shown  in  Figure  2. 

2.2  Material,  Geometry  and  Test  Methods 
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“13431  aV' 


Figure  1:  Load  History 


HIP  BETWEEN  CHANGE  IN  CRACK  LENGTH  AND  CRACK  LENGTH 


2.2.1  Material  and  Specimen  Geometry 

A  24  in.  (610  mm)  wide  center-cracked  panel,  manufactured  from 
0.182  in.  (4.6  mm)  thick  7075-T6  aluminum  alloy  was  used  (see  Figure  3). 
The  central  notch  was  1.663  in.  (42  mm)  in  length. 

2.2.2  Stress  Intensity  Factor  Coefficient 

Because  the  stress  intensity  factor  is  directly  related  to 
the  crack  growth  rate  Aa/AF,  the  f llght-by-flight  crack  growth  rate 
behavior  can  be  predicted  for  any  level  of  the  maximum  stress  intensity 
factor  •  This  functional  relationship  can  be  described  by: 

^  ~  ^^^max^  which  is  applicable  to  steady  state  spectra^^\  A  maxi¬ 
mum  stress  Intensity  factor  (K^  )  of  30  ksl  (33  MP  *^)  was  chosen 

to  correspond  with  the  maximum  applied  load  level  in  the  spectrum  and  was 
controlled  within  one  percent  through  load  shedding.  By  keeping  the  stress 
intensity  factor  "quasi"  constant,  a  determination  can  be  made  of  which 
distribution  function  Aa/AF  follows  at  a  particular  stress  intensity 
factor  level. 

2.2.3  Load  Shedding 

The  maximum  stress  intensity  factor  was  kept  "constant"  within 
one  percent  by  reducing  the  load  levels  for  every  0.050  in.  (1.27  mm) 
of  total  (2a)  crack  length.  This  was  determined  as  follows:  to  keep 
^a+Aa  percent  of  K^,  a  first  order  approximation  was  used. 


a+Aa 


^  1.01  K 

a 


Kg  +  ^  [(a+Aa) -a]  ^  1.01 

Aa  =  0.01  K 

a 


6 


da 


For  a  center-cracked  panel, 

P  /rra  /  . 

max _ aJ  secTra 

B  V  W 

w 


where  P  =  load;  B  ■■  thickness  of  plate;  and  W  ■»  width  of  plate.  A 

value  of  0.025  in  (.635  imn)  for  Aa  was  selected  as  the  growth  increment 

for  load  shedding  (Figure  2) .  For  most  of  -the  data  collected ,  the 

change  in  K  was  less  than  one  percent, 
max  ^ 

2.2.4  Test  Equipment 

A  500  kip(2.2MN)  static,  250  kip  (I.IMN)  dynamic  capacity 
load  frame  was  used  to  subject  test  specimens  to  flight-by-flight 
loading  previously  described.  The  load  frame  employed 
the  principles  of  closed-loop  servocontrol.  The  load  levels  and  cycle 
shape  were  stored  in  a  4096  byte  memory  digital  programmer  and  then  fed 
to  the  servo-  load  controllers.  The  500  kip  load  cell  connected  in 
series  with  the  test  specimen  provided  a  feedback  signal  for  load 
control. 

The  applied  test  loads  were  monitored  through  an  independent  data 
system  and  were  maintained  to  within  one  percent  of  the  programmed 
value.  The  crack  length  was  measured  using  a  binocular  zoom  microscope 
with  a  maximum  magnification  of  40X.  A  Mylar  scale  calibrated  in 
0.005  in  (0.127  mm)  divisions  was  attached  to  the  specimen.  Crack  length 
readings  were  estimated  to  be  within  +-002  inches  (.0508  mm).  This 
value  will  be  referred  to. as  the  measurement  precision,  e. 


K 

max 


a 


SECTION  3 


DISCUSSION  AND  RESULTS 

In  the  initial  test,  the  crack  length  was  measured  after  every 
fifteen  flights.  The  crack  length  versus  flight  data  was  tabulated 
and  put  into  a  statistical  analysis  program  provided  by  Purdue  University 
under  an  AFOSR  grant From  the  statistical  analysis  of  Aa/AF  a 
determination  was  made  as  to  which  of  four  (4)  statistical  distributions 
best  described  the  data.  The  data  are  shown  in  Figure  4. 

Crack  growth  rates  were  generated  using  the  secant  method  of 
differentiation.  The  differentiated  data  are  shown  in  Figure  5.  The 
resulting  Aa/AF  data  were  fit  to  four  distributions:  normal,  2  parameter 
lognormal,  3  parameter  lognormal  and  3  parameter  Weibull.  Of  the  four 
distributions,  the  3  parameter  lognormal  was  found  to  have  the  highest 
correlation  coefficient  for  all  the  data  sets  (Table  I  and  Figure  6). 
Because  of  the  high  correlation  coefficients  associated  with  the  four 
distributions,  all  appear  adequate  for  comparing  the  variability 
associated  with  one  set  of  crack  growth  rate  data. 

The  Aa/AF  data  from  the  first  panel  were  grouped  into  0.5  inch 
intervals  (Figure  7).  The  sample  means  of  the  groups  remained  fairly 
constant  as  the  crack  progressed.  The  means  were  indistinguishable 
from  one  another  using  the  student’s  t  test  at  a  954  confidence  level. 
These  four  groups  consist  of  all  the  data  collected  in  the  first  panel. 
Each  group  contains  between  forty  and  sixty  Aa/AF  observations.  Because 
the  sample  means  are  indistinguishable  as  the  crack  advanced,  a  distri¬ 
bution  can  be  found  which  describes  the  variation  Aa/AF  data  for  all 
crack  lengths  measured  for  a  given  number  of  flights. 
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Figure  5:  Constant  Control  Conditions.  Crack  Growth 
hate  Data  as  a  Function  of  Crack  Length  using  Secant 
Method  of  Differentiation 
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CUhULfiTIVE  STRNDflRD  NORMRL  PROBRBILITY  (PER  CEN 


3-PflRflMETER 

LOG  NORMAL  DISTRIBUTION  PLOT 


Figure  3-Parameter  Lognormal  Distribution 
Plot  ofaa/^^F  Data' 


96  %  CONFIDENCE  INTERVALS 


CRACK  LENGTH  A  (INCHES) 


The  largest  set  of  contiguous  points  was  used  to  evaluate  the  Aa/AF 
distribution.  The  sample  means  and  standard  deviations  are  shown  in 
Table  II.  Successive  points  were  then  removed  to  evaluate  the  effect 
of  increasing  the  crack  length  Increment  on  the  mean  and  variance. 

Three  differentiation  schemes  were  employed  to  investigate  the 
variability  in  the  Aa/AF  data.  The  sample  means  and  standard  deviations 
resulting  from  the  secant  and  two  seven  point  removable  strip  differen¬ 
tiation  methods  are  shown  in  Table  II.  When  "a"  versus  "F"  data  were 
fitted  with  a  polynomial  and  then  differentiated,  the  mean  remained 
constant;  the  variance,  however,  was  greatly  reduced  in  all  cases 
because  of  artificial  smoothing  Induced  through  this  process. 

By  removing  selected  points  and  differentiating  using  the  secant 

method,  the  effect  of  taking  larger  crack  growth  length  increment 

measurements  was  investigated.  Through  this  process,  the  sample  mean  of 

(9) 

Aa/AF  is  expected  to  remain  the  same.  Neulieb'  showed  that  for 
constant  stress  intensity  factor,  steady  state  conditions,  and  assuming 
independence,  the  coefficient  of  variation  (v)  should  decrease  as  one 
over  the  square  root  of  the  number  of  increments  which  were  added.  In 
this  study,  the  standard  deviations  as  shown  in  Figure  8  decrease 
faster  than  would  occur  for  independent  crack  growth  intervals.  It  is, 
therefore,  concluded  that  crack  growth  in  a  given  increment  may  be 
dependent  on  what  has  occurred  in  a  previous  increment.  An  interval 
of  fast  crack  growth  always  followed  by  an  interval  of  slow  crack  growth 
or  vice  versa  is  an  example  of  dependence.  In  the  case  of  dependent 
data,  the  coefficient  of  variation  will  decrease  at  a  faster  rate  than 
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Figure  8:  Comparison  of  TIeans  &  Std.  Dev.  for 
Successively  removed  Points 


than  the  square  root  of  N.  Dependency  may  be  attributed  to  actual 
material  behavior,  measurement  discrepancy  or  variability  in  load 
application. 

The  coefficients  of  variation  shown  in  Table  II  reach  an  acceptable 
level  of  under  twenty  percent  when  the  ratio  of  the  crack  length 
measurement  increment,  Aa,  to  the  measurement  precision,  e,  is  around 
ten  to  one.  A  Aa/e  ratio  of  10/1  should  be  considered  when  making  a 
determination  of  crack  growth  rates.  The  polynomial  methods  of 
differentiation  did  not  aid  in  the  assessment  of  variability  in  crack 
growth  rate  behavior. 

Determination  of  the  factors  effecting  the  variance  and  their 
relative  influence  would  require  further  study.  It  was  discovered 
that  a  different  value  for  the  variance  was  found  depending  on  which 
point  one  starts  with  in  the  differentiation  scheme  (Table  III).  This 
leads  to  a  suggestion  that  the  data  may  be  periodic.  This  periodic 
variability  is  another  indication  of  the  possible  dependency  of  the  data. 
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SECTION  4 


CONCLUSIONS 

It  was  found  that  the  flight-by-flight  load  history  used  for  this 
study  exhibited  steady-state  behavior.  This  was  determined  by  the 
smoothness  of  the  a  vs.  F  data  curve  and  the  repeatability  of  the  sample 
mean  and  standard  deviation  of  the  Aa/AF  data  as  the  crack  length 
increased.  Because  of  this,  a  distribution  function  can  be  found  which 
describes  the  Aa/AF  data. 

In  this  flight-by-flight  study,  the  three-parameter  lognormal 
distribution  function  with  a  correlation  coefficient  of  0.999  best 
described  the  Aa/AF  data.  The  two-parameter  lognormal  distribution 
function  with  a  correlation  coefficient  of  0.998  was  also  acceptable 
as  was  the  normal  distribution  function  with  a  correlation  coefficient 
of  0.991.  The  normal  distribution  function,  however,  allows  for 
negative  crack  growth  which  is  physically  unrealistic. 

The  method  of  differentiation  had  no  effect  on  the  mean  crack 
growth  rate.  The  variance,  however,  was  affected  by  the  method  of 
differentiation.  Using  a  polynomial  method  had  the  effect  of  decreasing 
the  variance  in  comparison  to  the  variance  attained  using  the  secant 
method  of  differentiation.  This  is  to  be  expected  because  the 
polynomial  method  is  a  smoothing  process.  The  secant  method  appears 
to  be  a  more  realistic  choice  in  establishing  crack  growth  rates. 

Increasing  the  crack  length  measurement  increment  size  had  no 
significant  effect  on  the  sample  mean.  The  variance  decreased  faster 
than  would  be  expected  for  independent  data.  Something  in  the  crack 
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growth  testing  process  gives  rise  to  dependency  in  the  data.  Possible 
contributors  are  the  measurement  technique,  load  shedding  process  and 
actual  material  behavior. 

The  coefficient  of  variation  is  related  to  the  relative  size  of 
the  measurement  precision.  A  ratio  of  crack  increment  size  to 
precision  (Aa/AF)  of  ten  to  one  gave  an  acceptable  coefficient  of 
variation  of  less  than  twenty  percent.  A  smaller  increment  size  gave 
a  higher  coefficient  of  variation.  Larger  increment  sizes  did  not 
improve  the  variability.  Increment  measurement  size  can  be  selected 
to  minimize  the  variability  in  f llght-by-flight  crack  growth  rate 
testing  by  ensuring  that  the  crack  length  measurement  interval  is 
large  compared  to  the  measurement  precision. 
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